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ABSTRACT 

An accurate picture of how free electrons are distributed throughout the Milky Way 
leads to more reliable distances for pulsars, and more accurate maps of the magnetic 
field distribution in the Milky Way. In this paper we test 8 models of the free electron 
distribution in the Milky Way that have been published previously, and we introduce 4 
additional models that explore the parameter space of possible models further. These 
new models consist of a simple exponential thick disk model, and updated versions of 
the models by Taylor & Cordes and Cordes & Lazio with more extended thick disks. 
The final model we introduce uses the observed Ha intensity as a proxy for the total 
electron column density, also known as the dispersion measure (DM). Since accurate 
maps of Ha intensity are now available, this final model can in theory outperform the 
other models. We use the latest available data sets of pulsars with accurate distances 
(through parallax measurements or association with globular clusters) to optimise the 
parameters in these models. In the process of fitting a new scale height for the thick 
disk in the model by Cordes & Lazio we discuss why this thick disk cannot be replaced 
by the thick disk that Gaensler et al. advocated in a recent paper. In the second part 
of our paper we test how well the different models can predict the DMs of these 
pulsars at known distances. We base our test on the ratios between the modelled 
and observed DMs, rather than on absolute deviations, and we identify systematic 
deviations between the modelled and observed DMs for the different models. For 
almost all models the ratio between the predicted and the observed DM cannot be 
described very well by a Gaussian distribution. We therefore calculate the deviations 
N between the modelled and observed DMs instead, and compare the cumulative 
distributions of N for the different models. Almost all models perform well, in that 
they predict DMs within a factor of 1.5-2 of the observed DMs for about 75% of 
the lines of sight. This is somewhat surprising since the models we tested range from 
very simple models that only contain a single exponential thick disk to very complex 
models like the model by Cordes & Lazio. We show that the model by Taylor & Cordes 
that we updated with a more extended thick disk consistently performs better than 
the other models we tested. Finally, we analyse which sightlines have DMs that prove 
difficult to predict by most models, which indicates the presence of local features in 
the ISM between us and the pulsar. 

Key v^rords: Galaxy: general; Galaxy: structure; ISM: structure; (stars:) pulsars: gen- 
eral 



1 INTRODUCTION 

Knowledge of the distribution of the free electrons in the 
Galaxy, rie, is important for several reasons. First, the de- 
tailed structure of rie is an essential ingredient for a full 
understanding of the multi-phase character of the Galactic 
interstellar medium (ISM). Furthermore, a model for the 
distribution of the thermal electrons is a crucial tool for 



estimating the distances of pulsars from their observed dis- 
persion measure 

(DMfl Finally, thr ough its integral DM, 
an Wo-model is required for an analysis of the strength and 
structure of the Galactic magnetic field from rotation mea- 
sures (RM), because RM is the integral over the line of 
sight of the product of the l.o.s.-co mponcnt of the magnetic 
field and rie (see, e.g., iHan et all i2006. : .Brown et al.ll2007l : 



* E-mail: Schnitzeler@mpifr-bonn.mpg.de 



^ The dispersion measure is the integral of the free electron den- 
sity along the line of sight, and is expressed in units of cm~^pc. 
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Figure 1. Galactic distribution of the isolated pulsars from Table [Tl (shown as blue circles) and pulsars in globular clusters from Table[2] 
(red circles). Filled circles indicate the 45 lines of sight that we selected in i; l4.2l the open circles indicate the remaining 23 lines of sight. 
In the left panel we show the sight directions in a Hammer- Aitoff equal area projection, the hatched band indicates the region |6| < 5°. 
The middle panel shows the positions of the pulsars and globular clusters projected onto the Galactic mid-plane; the background image 
shows on a linear grey scale the column densities from the NE2001 model that were integrated perpendicular to the mid-plane. In the 
right panel we show the distribution of the pulsars in the Rq, 2-plane, i.e. in the meridional planes through the Sun. 



Men et al.ll20()i:lNoutsos et al.ll2008l : IRov et al.l[200a : IValled 
20081 : iNota fc Katgertll20ld ). 

In our paper we analyse 8 models of the Galactic free 
electron density rie that have appeared in the literature. 
These m odels ranRe from s i ngle e x ponential thick disks, 
see, e .g., iBerkhuiisen et akl ( 20061 ). iBerkhuiisen fc Miilleil 
(|2008l ) and iGaensler et al.l (|2008l ). models that consist of 
a thin and a thick disk with di fferent radial and vertical 
scale heights (|G6mez et al.l 200ll). and the com plex, multi- 
com ponent models by iTavlor fc Corded (|l993l ). or 'TC93', 
and iGordes fc Lazid (|2002l ). also known as NE2001, that 
consist of an axisymmetric thin and thick disk, spiral arms, 
and a contribution by the Gum nebula (and, in the case 
of NE2001 additional components). We add four models to 
this list. The first model uses an exponential thick disk that 
is fitted to the latest available data on pulsars with accu- 
rate distances. We used a similar approach to fit new scale 
heights of the thick disk component in TG93 and NE2001, 
whilst maintaining the contributions by the other compo- 
nents in these models. With our final model we consider 
descriptions in which the prediction of DM is based on an- 
other observable quantity, viz. emission measure (EMfl. We 
then test the predictive quality of the available models using 
a new formalism that we developed. 

We arranged this paper as follows. We present and dis- 
cuss the data in Sect.[2l and in Sect.[3]we give an overview of 
the different models that we test. In Sect. |3] we describe the 
practicalities of modelling DM, in particular of the four mod- 
els we propose, and our criteria for selecting usable sight- 
lines. We then test the predictive quality of the available 
models in Sect. (5] Finally, in Sect. [6] we discuss sightlines 
for which the DMs prove very difficult to predict. Through- 
out this paper we have calculated statistics in a way that is 
robust against outliers. 



2 THE DATA 

Our analysis is based on a sample of pulsars with relatively 
accurate distances, either from a parallax measurement, or 
from an association with a Galactic globular cluster. In our 
analysis we selected lines of sight for which the fractional 
error in the parallax (isolated pulsars) or distance (pulsars 
in globular clusters) is less than 0.25; a s imila r approach was 
adopt ed bv lBerkhuiisen fc Miillerl (|2008l ) and lGaensler et al] 
(l2008h . Small distance errors make it easier to identify devi- 
ations between the modelled and observed DMs, which are 
the result of structure in the ISM that is not included in the 
model. 

Our sample consists of the (isolated) pulsars with 
known parallaxes that are listed on Shami Chatterjee's web- 
sit^fl (accessed on 27/7/2012), for which we took the posi- 
tions and me asured DMs from the A TNF pulsar catalogue 
(version 1.44; [Manchester et al ]|2005l . and thi s websitfEl) . For 
all the se pulsars we used the parallaxes that IVerbiest et al.l 
(|2010l ) corrected for Lutz-Kelker bias. There are no DMs for 

3 of the pulsars listed in the ATNF p ulsar catalogue (viz. , 
J0633-H746, J0720-3125, J1856-3754). iFreire fc We^ (|2010l ) 
published a parallax measurement for J 1738-1-0333 that was 
not yet in the list on Chatterjee's website when we accessed 
it, but we included it in our analysis. We corrected the paral- 
lax of J1738-I-0333 for Lutz-Kelker bias using the online tool 
at |http: / /psrpop.phys. wvu.edu/LKbias/ . This selection re- 
sults in a list of 41 pulsars with known DMs and parallaxes, 
where we counted the double pulsar J0737-3039 only once. 
The relevant properties of these pulsars are listed in Table 
[1] In the case of asymmetric parallax errors we made the er- 
rors symmetric by always using the highest of the two error 
estimates. For the 68 sightlines which we use in our analy- 
sis the difference between the upper and lower error limits 
amounts to at most 4%, and in the majority of cases to 
< 1%, of the bias-corrected parallax. Also, the percentual 
error in the DMs of these sightlines is at most 2%, which 



^ The emission measure is defined as the line-of-sight integral of 
the square of the free electron density, and is expressed in units 
of cm~®pc. 



http:/ /astrosun2. astro. corncll.edu/rcsearch/parallax/ 
* http:/ /www. atnf.csiro.au/research/pulsar/psrcat/ 
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means that the error in DM is smaller tha n the error in 
the pa rallax by about an order of magnitude. IVerbiest et al.l 
(j2010 l^ list an additional 14 pulsars with fractional parallax 
errors larger than 0.25 (six of these are larger than one). We 
did not include these 14 pulsars since adding them would 
considerably reduce the overall accuracy of the new sample 
compared to the old sample, and would increase the number 
of sightlines by only a modest amount. 

While the distribution of the parallax errors can be well 
approximated by a Gaussian distribution, the distribution of 
the distance — 1/ parallax is not Gaussian. It can be shown 
that distance follows the probability density function 



p(dist; TTobs, err^) 



1 



/2PT err^ dist^ 



where TTobs and err-^ are the mean and error of the parallax 
distribution. To avoid confusion with tt — parallax, we used 
'Pr to indicate the mathematical constant tt. In particu- 
lar, the most probable distance is not l/jTobs, but l/(0.5Trobs 
-I- 0.5-\/ TT^jjg -I- 8err|). The most probable distance is about 
11% smaller than l/7Tobs for the smallest TTobs/err^ = 4 
from our sample, and it approaches 1/TTobs asymptotically 
for large TToba/err^. The distance errors in Table. [T] define 
the most compact 68% confidence interval around the most 
probable distance. By using Lagrange multipliers it can be 
shown that this interval is bounded by distances a and b 
for which the probability density function at a equals the 
probability density function at b. 

Paulo Freire lists pulsars in globular clusters on his web- 
sit43 (accessed on 27/7/2012). We o btained the p roperties 
of their parent globular clusters from iHarrisI (| 19961 . and this 
websit^fl, version Dec. 2010). We list the globular clusters 
with their measured DM in table [2l When there is more 
than one pulsar in a globular cluster we calculated the aver- 
age DM and its 1-a spread. Clusters that are rich in pulsars, 
like NGC 104 and Terzan 5, show a statistically significant 
spread in pulsar DMs, which means that there is warm ion- 
ized gas inside globular clusters. The spread in DMs is how- 
ever small, only a few percent of the average DM for each 
cluster, so that we can average the pulsar DMs in each clus- 
ter, and use this DM as the characteristic value for that line 
of sight. This yields DMs for 27 globular clu sters with known 
distances. Following ICaensler et al.l (|2008l ) (who base their 
choice on a private communication of C. Heinke), we adopt a 
fractional uncertainty in the distance of each globular cluster 
of 15%^ 

ICaensler et al.l (|2008l ) include also pulsar DMs in the 
Magellanic Clouds in their analysis as upper limits on the 
Galactic DMs for these sightlines. These upper limits are 
much higher than the DMs we find for objects far from the 
Galactic plane, which means that they must contain a sig- 
nificant contribution from the Clouds themselves. Therefore 
we did not include these data in our analysis. 

In Fig. [T] we show the distribution of the 68 lines of 
sight from Tables [T] and [5] in 3 different projections. This 
figure shows that the sightlines towards pulsars with paral- 
lax distances sample structure in the Galactic ISM within 
a few kpc from the Sun, and the sightlines towards pulsars 

^ http:/ /www. mpifr-bonn.mpg.de/staff/pfrcirc/GCpsr.txt 
^ http://physwww.mcmastcr.ca/'~harris/mwgc.dat 



in globular clusters sample longer sightlines that mostly lie 
further than 3°- 4° from the Galactic plane. In ij 14.21 we de- 
scribe our criteria for selecting 45 sightlines from this sample 
that we then subsequently use to determine the parameters 
for the Tio models. In Sect. [5] we analyse how well the differ- 
ent models predict DMs for the subset of 45 'good' lines of 
sight, and for the full sample of 68 sightlines. 



3 OVERVIEW OF ELECTRON DENSITY 
MODELS 

In this section we give an overview of the models that were 
published in the literature, supplemented with a description 
(1) of the models M1-M4 that we introduced. Table |3] summa- 
rizes the inner workings of the rio-models for which we will 
analyse the quality of the DM predictions. In our overview 
we start with the simplest imaginable models, consisting 
of a single exponential thick disk, and build t oward s the 
multi-componen t s mo dels by iTavlor fc Corded l|l993l ) and 
ICordes fc Laiol l|2002l) . 

3.1 Plane-parallel exponential models with a 
single scale height, model Ml 

iBerkhuiisen fc MiiUeil (|2008D and IGaensler et"al] (|2008D 

have investigated models with a plane-parallel Ue- 
distribution with a single scale height. I.e., 

na{x, y, z) = nc{z) = noe"'^''''' 



so that the dispersion measure 



DPI{L, b) 



Tic dl — 



noh 



[1 



-Lsm(t6|)/h, 



(2) 



(3) 



where L is the length of the line of sight, and b is its 
Galactic latitude. The fraction on the right hand side is the 
DM of an infinitely long line of sight, DMoo, and the term 
between brackets corrects DMoo for a finite path length. 
For a uniform free electron density in the Galactic mid- 
plane the DM that is accumulated perpendicular to the 
disk can be easily calculated as DMx = DMsin(|&|). In 
this paper we will often refer to this quantity, and to the 
DM_L of an infin i tely l o ng line of sight, DMq q i = noh. 
Berkhuiison et al' ('2006') . IBerkhuiisen fc Miillej (j2008l '). and 
Gaenslcr ot al. (2008) find somewhat different DMoo,± of 25, 
22, and 26 cm^'^pc resp., which can be related to the differ- 
ent criteria that these authors used to select lines of sight 
for their analysis (we list these in Table |3| . 

Our model Ml consists of only a single exponential 
thick disk, similar to the models by Berkhuijsen et al. and 
Gaensler et al., for which we determine the best scale height 
in § 14.31 based on the data from Tables [T] and (2] 



3.2 Plane-parallel 2-component models with a 
radial scale length 

I Gomez et al.l (|200ll ) considered models that consist of a thin 
and thick disk that each show also a dependence on Galac- 
tocentric distance. They determined the parameters of these 
models using a data set of 109 pulsars. (Throughout our pa- 
per, when we use the term 'Galactocentric radius' or 'Galac- 
tocentric distance' we are referring to a distance projected 



4 D.H.F.M. Schnitzeler 



Table 1. Properties of the 41 pulsars in our sample with fractional parallax errors < 0.25. For each pulsar we list its Equatorial and 
Galactic coordinates, its measured DM, and the distance from the Sun (with associated error) and height above the Galactic mid-plane 
that were derived from parallax measurements. *: these distance errors define the most compact 68% confidence interval around the most 
probable distance. We excluded pulsars from our analysis for the following reasons: °: pulsar lie s within 5° o f the Galactic plane; line 
of sight is affected by strong Ho emission in its vicinity; line of sight has a BRT_OBJ flag in iFinkbeinej ||2003| ). In S 14.21 we describe 
why we leave out sightlines that have an 'a','b' or 'c' flag, and wc use the remaining 27 sightlines for our analysis. 
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onto the Galactic mid-plane) For only 8 of their pulsars a 
parallax distance was available, and for 27 of their pulsars 
the distance is based on an association with either a globu- 
lar cluster or a supernova remnant. Most of the remaining 
distances are kinematic. The errors in the latter were dis- 
cussed by iGomej (|2006l ) and were found to be significant, 
especially at the positions of the spir al arms. The genera l 
form of the electron density from which [Comez et"aL (|200ll ) 
calculated DMs is 



ne(rGC, z) = ni 



f{rQ/hr,i 



/(rGc/fer,2) 

f{rQ/hr,2) 



z — Lsin(|fo|), and f{x) = e ^ (which we refer to in Ta- 
ble [3] as 'GBCa') or f{x) = sech2(x) (^'GBCb'). In Table B 

200 



we list the fit parameters obtained bv lGomez et al , . _ _ . 
Their model is somewhat reminiscent of that bv lLvne et al.l 
which consisted of a thick and thin disk component 
(the latter having a small vertical sca le height). The thick 
disk in the model bv lLvne et al.l (jToSS') provides a constant 
free electron density, contrary to the other models that we 
tested. 



(4) 



where z is the height above the Galactic mid-plane, 
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Table 2. Properties of the 27 globular clusters in our sample that are known to contain at least one pulsar. For each line of sight we list 
its Equatorial and Galactic coordinates, its measured DM, the number of pulsar the globular cluster contains, and its distance from the 
Sun and height above the Galactic mid-plane. We excluded lines of sight for the following reasons: line of sight lies wit hin 5° of th e 
Galactic plane; line of sight is affected by strong H« emission in its vicinity; line of sight has a BRT_OBJ flag in Finkb eTneJ l l2003f) . 
In ij. 14.21 we describe why we leave out sightlines that have an 'aVb' or 'c' flag, and we use the remaining 18 sightlines for our analysis. 



Cluster name RA2000 o DEC2000 ' b DM ^psi dist \z\ 

[°] [°] [cm-3pc] [pc] [pc] 



NGC 104 


00:24:05.2 


-72:04:51 


305.90 


-44.89 


24.37 


23 


4500 


3200 


NGC 1851 


05:14:06.8 


-40:02:48 


244.51 


-35.03 


52.15 


1 


12100 


6900 


NGC 5024 


13:12:55.3 


-1-18:10:06 


332.96 


79.76 


24.00 


1 


17900 


17600 


NGC 5272 


13:42:11.6 


-f 28:22:38 


42.22 


78.71 


26.37 


4 


10200 


10000 


NGC 5904 


15:18:33.2 


-1-02:04:52 


3.86 


46.80 


29.49 


5 


7500 


5500 


NGC 5986 


15:46:03.0 


-37:47:11 


337.02 


13.27 


92.17 


1 


10400 


2400 


NGC 6121'' 


16:23:35.2 


-26:31:33 


350.97 


15.97 


62.86 


1 


2200 


600 


NGC 6205 


16:41:41.2 


-1-36:27:36 


59.01 


40.91 


30.17 


5 


7100 


4700 


NGC 6266 


17:01:12.8 


-30:06:49 


353.57 


7.32 


114.08 


6 


6800 


900 


NGC 6342 


17:21:10.1 


-19:35:15 


4.90 


9.72 


71.00 


1 


8500 


1400 


NGC 6397 


17:40:42.1 


-53:40:28 


338.17 


-11.96 


71.80 


1 


2300 


500 


Terzan 5" 


17:48:04.8 


-24:46:45 


3.84 


1.69 


237.82 


34 


6900 


200 


NGC 6440'^ 


17:48:52.7 


-20:21:37 


7.73 


3.80 


222.80 


6 


8500 


600 


NGC 6441'''= 


17:50:13.1 


-37:03:05 


353.53 


-5.01 


232.24 


4 


11600 


1000 


NGC 6517 


18:01:50.5 


-08:57:32 


19.23 


6.76 


182.45 


4 


10600 


1200 


NGC 6522" 


18:03:34.1 


-30:02:02 


1.02 


-3.93 


193.03 


3 


7700 


500 


NGC 6539 


18:04:49.7 


-07:35:09 


20.80 


6.78 


186.38 


1 


7800 


900 


NGC 6544" 


18:07:20.6 


-24:59:50 


5.84 


-2.21 


134.00 


2 


3000 


100 


NGC 6624 


18:23:40.5 


-30:21:40 


2.79 


-7.91 


87.63 


6 


7900 


1100 


NGC 6626 


18:24:32.8 


-24:52:11 


7.80 


-5.58 


120.10 


12 


5500 


500 


NGC 6656 


18:36:23.9 


-23:54:17 


9.89 


-7.55 


91.40 


2 


3200 


400 


NGC 6749" 


19:05:15.3 


-f01:54:03 


36.20 


-2.21 


192.85 


2 


7900 


300 


NGC 6752 


19:10:52.1 


-59:59:04 


336.49 


-25.63 


33.36 


5 


4000 


1700 


NGC 6760" 


19:11:12.0 


-1-01:01:50 


36.11 


-3.92 


199.68 


2 


7400 


500 


NGC 6838" 


19:53:46.5 


-1-18:46:45 


56.75 


-4.56 


117.00 


1 


4000 


300 


NGC 7078 


21:29:58.3 


-1-12:10:01 


65.01 


-27.31 


66.88 


8 


10400 


4800 


NGC 7099 


21:40:22.1 


-23:10:48 


27.18 


-46.84 


25.08 


2 


8100 


5900 



Table 4. Parameters of the models GBCa (which uses f(x)=exp(- 
x)) a nd GBCb (f(x)=sec h^(x)) that give the best fit to the data 
from lGomez et al.l 1 I2OOII) . 



Model 


f(x) 


n(i'=r0,z=O) 




hr 






[cm^^] 


[kpc] 


[kpc] 


GBCa 


exp(-x) 


2.03 X 10-2 


1.07 


30.4 






0.71 X 10-2 


0.05 


1.5 


GBCb 


sech^ (x) 


1.77 X 10-2 


1.10 


15.4 






1.07 X 10-2 


0.04 


3.6 



3.3 



Multi-component models TC93 and 
NE2001a,b,c, models M2 and M3 



This family of models is built around the same basic 
structure for the Milky Way, which was first described by 
iTavlor fc Cordej (|l99lj ) . Their model consists of a thick disk 
with a large radial scale length and a large extent perpen- 
dicular to the mid-plane, an additional thin disk component 
close to the Galactic pla ne, 4 spiral arms whose locatio ns are 
based on the model by iGeorgehn fc GeorgelinI (| 19761 ) that 
were subsequently modified slightly, and an enhancement 
of the free electron den sity i n the Gum nebula. The mod- 
els bv lLvne etaP (jl985l ) and lCordes et all (|l99ll) contained 
only t hin d isk and thick disk compo nents, but Ghosh fc Rad 
l|l992l ') and lTavlor fc CordesI l|l993l ') argued the need for spi- 



ral arms. In the remainder we will refer to the model by 
iTavlor fc CordesI (|l993ll as 'TC93'. 

The model that was defined bv lCordes fc Laziol (|2002l ) 
consists of 7 components, each with its analytical expression 
for the free electron density rio as a function of position in 
the Galaxy. There are 3 smooth components: the spiral arms 
and the thin and thick disk. Together, these are defined by 
28 parameters. Then there is a Galactic-centre component 
(defined by 3 parameters), and a description of the Local 
ISM (which requires 36 parameters). Finally, there is a sig- 
nificant number of clumps and voids with which the authors 
describe the smaller-scale structure in the ISM. 

We compare the DM-predictions of two versions of the 
NE2001 model with the data in Tables [T] and [1 The first 
version is defined by the web interfac^ll (NE2001a); the sec- 
ond one is defined by the Fortran code that can be down- 
loaded from the same web address (NE2001b). In Sect. [5] 
we will show that the web interface and Fortran code pre- 
dict different DMs for the same sightline s. Three papers 
have recently been published (|Sun fc Reichll2oTol : iLaw et al] 
l201ll :|j ansson fc Farradlioil ) in which the authors replaced 
the scale height of the thick disk in the NE2001 model 



http: / /rsd- www. nrl.navy.mil/7213/lazio/ne_model/ 
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Table 3. Overview of the models that we considered, and the accuracy of the DMs that they predict. For each model we list the factor 
N by which the modelled DM deviates from the measured DM (see Sect. (5)1 for the 50th, 75th and 90th percentile of the ensemble of 45 
selected lines of sight, and all 68 lines of sight, as we describe in Sect.fS] We used a linear interpolation between data points to calculate 
the N that correspond to exactly the 50th, 75th and 90th percentiles. *: N was not calculated for model M4 when all 68 sightlines are 
use d, because the EM for som e of these sightlines cann ot be calculated reliably (see ^H!T|. References: (1) |Bcrkhuiiscn & Miillcr (200^; 
(2'llGaensler et al.l ||2008|'I: O'llGomez et al'l ll200ll 1: (41 lTavlor fc Corded l ll993h : (•5 HCordes &: Lazid l|2002l') : (61 .Sun fc Reich L2010) : (7) 
iLaw et al. I 1I2OIII I; (81 Ijansson fc FarraJ \201'£ ~^ 



Model 


Characteristics 






Predictive quality N for 










50% 


75% 


90% 


50% 


75% 


yu/o 








of 45 lines of sight 


of all lines of 


sight 


Exponential models with a single scale height 
















BM08'' 


h = 0.93 kpc, DM00, ± = 21.7 cm-^pc 




1.33 


1.48 


1.74 


1.41 


1.65 


2 28 


GMCM08'' 


h = 1.83 kpc, DMoo^± = 25.6 cm-^pc 


2 


1.20 


1.70 


2.12 


1.49 


2.09 


3.13 


Ml^^ 


h = 1.59 kpc, DM00, ± = 24.4 cm-^pc 


this paper 


1.21 


1.60 


2.03 


1.39 


1.94 


2.86 


Plane-parallel 2-component models with a radial scale length 
















GBCa 


Eqn.|4]with f{x) = e~^ 


3 


1.32 


1.46 


1.72 


1.36 


1.63 


2.17 


GBCb 


Eqn. Hlwith f{x) = sech^(x) 


3 


1.35 


1.48 


1.68 


1.37 


1.60 


2.15 


Multi- component models 
















TC93 


Model by Taylor & Cordes 


4 


1.33 


1.50 


1.81 


1.35 


1.56 


2.33 


M2<i 


TC93 with a thick disk of 1.59 kpc 


this paper 


1.20 


1.37 


1.82 


1.22 


1.53 


2.03 


NE2001a 


Model by Cordes & Lazio (Web interface) 


5 


1.26 


1.79 


2.29 


1.29 


1.86 


2.86 


NE2001b 


Model by Cordes & Lazio (Fortran routine) 


5 


1.25 


1.76 


2.25 


1.27 


1.77 


2.85 


NE2001C 


Model by Cordes & Lazio 


6,7,8 


1.25 


1.52 


1.85 


1.32 


1.62 


2.04 




(Fortran routine + thick disk from GMCM08) 


















NE2001b with a thick disk of 1.31 kpc. 


this paper 


1.29 


1.51 


1.83 


1.35 


1.59 


2.02 




DM00, x= 20.5 cm-3pc 
















Structured model using DM(EM) 
















M4d 


DM00 = 22.67 EM^64 for extrapolation 


this paper 


1.98 


2.53 


3.24 


* 







to finite distance 



^ : constrained by 38 lines of sight, errdist/dist < 0.5 , effective thq < 1 

^ : constrained by 16 lines of sight with |fe| > 40°, err^iat/dist < 0.33 

: constrained by 22 lines of sight from Tables [T] and [2] with |6| > 32.5°, err^iat/dist < 0.25, effective rnct < 1- 
: constrained by 44 lines of sight from Tables [l] and [2] with |fe| > 5°, errdist/dist < 0.25, effective t-hc, < 1. 



(which was h = 970 pc[f|) wit h the much lar g er va lue of 
h = 1830 pc that was found bv lGaensler et all (|2008| 1. The 
mid-plane dens i ty wa s also replaced by the value used by 
iGaensler et al.l (|2008h . which changes the DM00, ± of the 
thick disk from 33 cm~^pc down to 26 cm~^pc. We will 
refer to this as model NE2001c. 

With our models M2 and M3 we investigate whether 
changing the scale height of the thick disk in TC93 resp. 
NE2001b improves their accuracy in predicting DMs. (i) l4.4p 

3.4 M4: A structured model based on the 
DM(EM)-relation 

The strength of the models we discussed in the previous sub- 
sections is that they require only a small number of model 
components. However, on small scales the ISM is not as 
smooth as these models predict. With our model M4 we 
test whether one can predict DM more accurately by us- 
ing the emission measure EM as a proxy to predict DM00, 
and then use one of the other models to correct for the fi- 
nite distance to the pulsar. This approach could in principle 

^ The Fortran code uses 970 pc, but Table 3 in lCordes fc Laziol 
l|2002l l lists 950 pc. We will adopt 970 pc as the correct value. 



outperform other models, since the EM can reflect the struc- 
ture in DM on small physical scales. In § 14.51 we will test 
whether model M4 indeed gives a smaller scatter in DMoo^x 
than other models. 

One caveat in this approach is that the same electrons 
are as sumed to produce both DM and EM. ISnowden et al.l 
(1997') showed by analysing ROSAT X-ray data that the 
mid-plane density of free electrons in the Hot inter-cloud 
medium (HIM) is 3.5 x 10~^ cm~'^, which is an order of 
magnitude smaller than the mid-plane density of the free 
electrons in the Warm ionized medium (WIM). We there- 
fore neglected their contribution to DM, and, as a conse- 
quence, both the observed DM and EM are produced only 
by electrons in the WIM. 



4 DETAILS OF THE INVESTIGATED 
MODELS 

In this section we describe the practicalities of calculating 
the model DMs. We will first describe how we calculate EM, 
since we will use EM to remove sightlines in the direction 
of complex Ho? features. In § 14.21 we then describe how we 
selected 45 'clean' lines of sight from Tables [1] and [2] that 
we will use to fit model parameters. Finally, in § 14.31 to 
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14.51 we will discuss the inner workings of models M1-M4 
in more detail. In § 14.41 we will also explain why the thick 
disk from NE2001 cannot be exchanged for the thick disk 
from GMCM08. 

For the models from the literature (TC93, NE2001a,b,c, 
BM08, GMCM08, GBCa and GBCb), we took the param- 
eters such as vertical scale height h, mid-plane density no, 
and radial scale length (for GBCa and GBCb) from the re- 
spective papers, to predict DMs for each of the 45 lines of 
sight from Tables [T] and [21 



4.1 Calculating EM 

The EMs are calculated in two steps. In the first step, the 
observed Ha intensities are translated into 'observed' EM, or 
EMobs- In the second step these are converted to EM for the 
entire line of sight, by correcting for interstellar extinction. 

The first step in this calculation s tarts with the all- 
sky map produced by iFinkbeineJ (120031') from the WHAM , 
SHASSA and VTS S su rveys (see | Haffner et al] l|2003h . 
iGaustad et"all l|200ll ) and lDennisonet al. (199^) resp.). The 
requ ired conversio n facto r is quite well established; see, 
e.g., Haffner et al. l|l998h . and the detailed discussion 



IValls-GabaudI l| 19981 ) . in particular his 'Case B' scenario of 



optically thick Lyman continuum radiation. It appears that 
different authors agree on the conversion factor EMobs /^Ha 
to within 2% for temperatures around 8000 to 9000 K, 
which are usually assumed for this ISM phase l|Revnoldsl 
1 19851 : iMadsen et a l. 2006 ). We have therefore used the 
relation given by iHaffner et all (|l998l ). viz.: EMobs = 
2.75 (T/10'')°'^/hc<, with EM in the usual units, Ihc in 
Rayleigh, and T = 8000 K. 

The second step, i.e. the correction for extinction, is 
less direct. The correction is e'^, where r is the effective 
optical depth of the Hg-emission . The total op t ical d epth 
foUows from the E{B - V) that ISchlegel et al] (|l998l ) de- 
rived, assuming R — Ay / E{B — V) = 3 .1, and taking 
Ahc/Av = 0.81 (see lDickinson et al.|[20o3 . and references 
therein). The effective optical depth depends on the relative 
distribut ions of the diffuse gas and the dust along the line 
of sight; iDickinson et al.l (|2003h derived that the effective 
optical depth of the Ha line is 0.33 ± 0.10 times the total 
optical depth. The correction factor for interstellar extinc- 
tion then becom es e*"*'^^ ^ e(b-v) did no t use the values 
of Av given bv iRowles fc FroebrichI ||2009D . because their 
results are quit e comparable to those of Schlegel et al. for 
|6| > 5°, while IPeek fc Graved (|2010l ) confirm the Schlegel 
et al. E{B — l/)-values to within 0.05 magnitudes, although 
the maps that Peek et al. made have a resolution of only 
4.5°. 

We limit t he an alysis to pulsars with \b\ > 5°, because 
ISchlegel et"all (|l998 'l warn that the E{B - V) they derive 
are on ly re liable for these latitudes. A lso, IDickinson et al.l 
l|2003l ) and lBerkhuiisen fc Miilleil { 20081 ) point out that one 
can only reliably correct for interstellar extinction when the 
effective optical depth is smaller than 1, which excludes most 
lines of sight with < 5°. 

There is an additional complication when estimating 
EM for pulsars in globular clusters. Globular clusters are 
bright in Ha, but most of this emission is produced by the 
stars in the cluster, not by the gas. We can therefore not use 



the Ha intensity of the cluster itself to calculate the Galactic 
EM along the line of sight. Instead, we calculate the average 
EM for a 3x3 rectangular grid of lines of sight, excluding the 
central line of sight. The cells in this grid measure 18'xl8', 
which puts the outer lines of sight outside the edge of the 
largest globular clusters in our sample. 



4.2 Selection of the lines of sight 

We selected a subset of 27 lines of sight from Tableland 18 
lines of sight from Table [2] to determine the fit parameters 
for models M1-M4 from Table [3] First we left out all lines of 
sight that lie within 5° of the Galactic plane. This selection 
avoids regions close to the Galactic plane where structure in 
the ISM generally becomes very complicated, and difficult to 
describe without detailed modelling. Also, as we described 
in § 14.11 this avoids lines of sight with unreliable EM where 
our model M4 cannot be used. 

For the remaining lines of sight we added 2 more se- 
lection criteria. The high-resolution Ha maps from VTSS 
and SHASSA showed strong, structured Ha emission in 
the direction of the pulsars J0659-I-1414, J1643-1224, and 
J2055-I-3630, and the 2 globular clusters NGC 6121 and 
NGC 6441. The lines of sight towards J1713-(-0747, J1909- 
3744, and NGC 6441 furthermore have 'BRT _ OBJ' fiags in 
the Ha mask that was created bv lFinkbeineil l|2003l ). which 
indicates regions where the automated removal of bright 
point sources cannot be trusted. Since both effects pro- 
duce unreliable EM we excluded these lines from our sample 
(NGC 6441 is removed for both these reasons). Finally, we 
checked that the remaining lines of sight all have an effective 
optical depth < 1, which is required to reliably correct for 
interstellar extinction, as described in § 14.11 

The remaining sample of lines of sight that we will use 
to determine the parameters of the models M1-M4 contains 
27 lines of sight towards isolated pulsars, and 18 lines of sight 
towards pulsars in globular clusters. As a consequence of our 
choice to use relatively unsophisticated selection criteria, a 
few of the selected lines of sight will have DMs that are 
difficult to model, as will become clear in the remainder of 
this paper. We will return to this in Sects. [5] and [6] 



4.3 Parameters for model Ml 

As described in Sect. [21 distances derived from parallaxes 
follow a complicated, non-Gaussian error distribution, which 
needs to be taken into account properly when fitting a scale 
height to these data points. We find the best-fitting scale 
height by maximising the likelihood C, where each pulsar 
with a parallax measurement contributes a probability den- 
sity for the distances given by Eqn. [U and each globular 
cluster with a DM measurement contributes a Gaussian dis- 
tribution of the distances. The log(likelihood) is given by 



log(£) « ^ -i 



L_prcd(h) 



pulsars 



err_7r^ 



2 log(L_pred(h)) 



1 {L - L_pred(h))^ 

2 err_L2 



(5) 



clusters 

where tt and L are used as shorthand for the parallax and 
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1000 r 



100 r 




DM sln(lbl) [crrr^pc] 



Figure 2. Distribution of height above the Galactic mid-plane 
\z\ vs. DMx = -DM sin{|b|) for selected pulsars from Table [Tl and 
globular clusters from Table [2] (crosses and filled circles resp.). 
The dotted line indicates the asymptotic value for DMsin(|6|) = 
24.4 cm~^pc that we calculated from all lines of sight with \z\ > 
4 kpc, which are indicated with black symbols. Colours iden- 
tify lines of sight that belong to the same latitude bin, and the 
coloured lines represent the best fits to the points with the same 
colour. The 5 colours shown in this Figure correspond to the data 
points that are drawn with solid error bars in Fig. |3] red belongs 
to the latitude bin 5° < |6| < 10°, dark blue to the next latitude 
bin, etc. We indicate the pulsar J0630-2834 that we did not use 
to determine h, but that we did use in our analysis from Sect. \5\ 



¥- ¥=- 



10 20 30 40 50 60 70 

abs(Gal. latitude) [°] 

Figure 3. The distribution of fitted scale heights h as a function 
of Galactic latitude. The vertical error bars reflect the 1-cr spread 
in the derived h, and the horizontal error bars indicate the lati- 
tude range Ab of the lines of sight that are used in the fit. The 
rightmost data point covers all latitudes |6| > 50°, and is centred 
on the average b of this latitude range, weighted by the surface 
area of each annulus db. With this choice of bin widths there are 
at least 5 lines of sight per latitude bin. Points with dashed lines 
use the same Ab as the data points immediately to the left of it. 



distance of the pulsar/globular cluster resp., I/_pred(h) in- 
dicates the distance of the pulsar/globular cluster that is 
predicted from its observed DM and input scale height h, 
and 'err' stands for the measurement error. We numerically 
evaluated Eqn. [5]for a grid of scale heights, and we fitted a 
parabola to the maximum of the distribution of log(yC) as a 
function of scale height to find the overall best-fitting scale 
height. We used this technique to determine the optimum 
scale height for all three models M1-M3 that we introduced. 

GCMC08 previously showed that the value of the scale 
height depends on which sightlines are used to fit the expo- 
nential thick disk to. While they changed the minimum |6| 
to find the best-fitting scale height for the thick disk, we di- 
vide the sightlines over bins in Galactic latitude and analyse 
how the fitted scale heights differ between the latitude bins, 
see Figs. [2] and [31 

First, when fitting a scale height, we solve for the height 
above the Galactic midplane 2 as a function of DMx in- 
stead of the other way round, since the errors in z are larger 
than the erro rs in DMx by about an order of magnitude 
(or more). As Savage, pointed out this can 

lead to a biased estimate of the fitted scale height, and we 
will return to this in the final paragraph of this section. The 
curves that we fitted to the data points are required to con- 
verge asymptotically to the DMx at infinity, DM00, ±, in all 
latitude bins. We approximated DM00, ± by the DMx of the 
8 lines of sight with |z| > 4 kpc, on the assumption that very 
little DMx is built up beyond those \z\. These 8 data points 
consist of 7 globular clusters and 1 pulsar that all lie at least 
27° from the Galactic plane, and are indicated in black in 
Fig. [1 This way we find DM^o.x = 24.4 ± 4.2 cm'^'pc. We 
then determined the scale height of the exponential thick 
disk for each latitude bin by maximising log(£) in the way 
we described in the first paragraph of this section. 

We colour-coded sightlines that belong to the same lat- 
itude bin, and the colours in Fig. [5] match the colours in 
Fig. [3l We used a Monte Carlo simulation with 10* draws 
to simulate the effect that the scatter in DM00, ± has on the 
fitted scale heights, and we will use the same approach to 
determine the uncertainty in the fitted scale heights in mod- 
els M1-M3. The horizontal error bars in Fig. [3] indicate the 
range that is covered by each bin in Galactic latitude, the 
widths of these bins were chosen in such a way that each bin 
contains at least 5 sightlines that we can fit a scale height 
with. Data points with dashed error bars in Fig. [3] show the 
fit results when the bins are shifted further off the Galactic 
plane by half a binwidth; adding these bins better samples 
the fitted scale height /i as a function of Galactic latitude. 

As Fig.|3]shows, the scale heights that are fitted to data 
points in the higher latitude bins are consistent with each 
other, only the scale height fitted to the blue data points 
in Fig. [5] (sightlines with 10°< |fe| < 25°) jumps out. When 
fitting the scale height of the exponential thick disk in our 
model Ml we stayed clear of this region by conservatively 
selecting only sightlines with |fo| > 32.5°; this selection also 
avoids sightlines between 17.5° < |6| < 32.5° where the tran- 
sition to the smaller scale height that was fitted to the blue 
data points occurs. We found a scale height h = 1586 ± 295 
pc f or the thick disk in Ml. 

ISavage fc Wakkeij (|2009l ) point out in the Appendix of 
their paper that an exponential model cannot reproduce ob- 
served DMx that are larger than the DM00,x of the model. 
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Figure 4. Comparison between the thick disk components in models M1-M3 and the data these components were fitted to. Crosses 
indicate isolated pulsars, filled circles indicate pulsars in globular clusters. Black data points show the original DMsin(|6|) and \z\ data 
points from Tables [T] and (2] and the black model curve shows the exponential disk model Ml that was fit to these data. Red data points 
show the DMsin(|6|) residuals (explained in ij 14.41 1 and model curves for the exponential disk components in models M2 and M3 that 
were fitted to the red data points. The model curves asymptotically approach the vertical dotted line, which is determined from the 
DMsin(161) (model Ml) or residual DMsin(|f)|) (M2,M3) of target sources at \z\ > 4 kpc. 



and that this introduces a bias in the derived scale height of 
the model. The y conclude that this ex plains why the scale 
height found by lGaensler et aTl l|2008l ) is much larger than 
the 'canonical' value that was found in t he literature up till 
then. Although the selection effect that I Savage fc Wakken 
l|2009l ) describe is real, Fig. [2] shows that it affects only 
a few of the 45 lines of sight that we used in our Monte 
Carlo simulation: only the sightlines towards NGC 5986, 
NGC 6517, and NGC 6539 have DM^ that lie within 1 a 
(=4.2 cm~^pc) of DMoo,x=24.4 cm~^pc. Smaller asymp- 
totic values for DMoo,x are not often encountered in our 
Monte Carlo simulation. Furthermore, lines of sight far from 
the Galactic plane (]6| > 40°) consistently reproduce scale 
heights of about 1.6 ± 0.3 kpc in our Monte Carlo simula- 
tion, whi ch is comparable to the scale height h = 1.4llQ'2f 
kpc that ISavage fc Wakkeij (|2009l ') themselves find for these 
latitudes. We therefore conclude that the scale heights we 
derive by fitting DMx as a function of 1^1 instead of the 
othe r way round are no t stron gly affected by the fitting bias 
that lSavage fc WakkeJ (|2009l ') identified. 

4.4 Parameters for models M2 and M3, 
discussion of NE2001c 

As Figs. [2] and [3] showed, an exponential disk model with 
a single scale height cannot provide the best fit at all lati- 
tudes, and the blue and red data points in Fig. [2] that in- 
dicate sightlines within 25° of the Galactic plane show con- 
siderable scatter. Furthermore, regions with extended Ha 
emission like the Gum nebula are known to add significant 
DM to pulsars that lie beyond it. Summarising, the Galac- 
tic ISM contains components in addition to the exponential 
thick disk that we considered so far. In this section and the 
next we investigate ISM models that model some of these 
additional components. Since we will follow the same ap- 
proach to construct models M2 and MS we will discuss both 
models here. 

Since the exponential thick disk in model Ml has a 
much larger scale height than the thick disk in TC93 and 
NE2001 (1.6 kpc vs. 0.9 kpc) we will investigate whether 
changing the thick disk in TC93 and NE2001 can improve 



the model performance in predicting DM. To remain inter- 
nally consistent with the other components in these models 
we first subtract the DM contributions by all components 
in TC93 and NE2001 except their thick disks (which we 
will refer to as the 'DM residuals') from the observed DMs. 
We then fit a new exponential thick disk model to the DM 
residuals using the procedure described in § 14.31 Contrary 
to Ml, where we used only sightlines at |&| > 32.5°, we use 
all sightlines at |ti| > 5° to determine the scale height of 
the thick disk in models M2 and MS: we assume here that 
TC9S and NE2001 model the additional ISM components 
perfectly, so that subtracting their predicted DMs from the 
observed DMs (= our definition of the residual DM) leaves 
only the contribution by the exponential thick disk. This 
solution is however not perfect. Since the relative magni- 
tudes of the DM contributions by the thin and thick disk 
(for example) follows from a complex optimisation proce- 
dure, changing the contribution of one without changing the 
contribution of the other can still lead to some internal in- 
consistencies. 

In Fig. 0] we compare the z vs. DM^ predicted by MI- 
MS using the best fits to the data points (Ml: black data 
points) and DM residuals (M2, M3: red data points). The 
black line indicates the best model fit of Ml in all S panels, 
and the red lines show the best fit to the data by model 
M2 (middle panel) and MS (right-hand panel). The verti- 
cal dotted line indicates the asymptotic value of DMoo.x, 
which is the same for M2 as for Ml. The reason for this is 
that the thin disk, Gum nebula and spiral arms in TC93 
do not contribute to the DMs of the 8 sightlines that we 
used to calculate DMoo,x. The local ISM components and 
local spiral arm that are included in NE2001 do contribute 
to the DMs of these 8 sightlines, so that correcting for their 
DM contributions leads to a smaller DMoo.x = 20.52 ± 2.72 
cm^'^pc. 

The scale height that we found for the thick disk in M2 
is 1587 ± 309 pc. Since the local ISM components and local 
spiral arm in NE2001 contribute to the DMs of sightlines at 
|6| > 40°, while GMCM08 assumed that only the exponen- 
tial thick disk component contributes DM to sightlines this 
far from the Galactic plane, the exponential thick disk from 
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Figure 5. The DM at infinity predicted by the exponential thick 
disk in model TC93 for 6=0° , and the prediction by the exponen- 
tial thick disk in model M2 when a radial cut-ofT at 27.7 kpc is 
used. We determined this cut-ofT radius by minimising the vertical 
separation between the two model curves. 

GMCM08 cannot replace the original thick disk in NE2001. 
The exponential thick disk that we fitted to the DM resid- 
uals for NE2001 has a scale height of 1313 ± 181 pc, which 
is in fact smaller than the scale heights that we found for 
Ml and M2, and much smaller than the scale height found 
by GMCM08. 

To find a physically reasonable radius for the exponen- 
tial thick disk that we fitted in M2, we require that this thick 
disk predicts the same DM at infinity as the thick disk in 
TC93. (We only do this for model M2, since we will show in 
Sect. [5] that M2 predicts DM best of all the models that we 
tested) We first determined the DM predicted at infinity by 
the thick disk in TC93 for a range of Galactic longitudes, 
at 6=0° . We then minimised the mean square difference be- 
tween the DM at infinity predicted by the thick disk com- 
ponents in TC93 and M2 for a range of trial radii. By fitting 
a parabola to this minimum we found that a radius of 27.7 
kpc gives the overall best match between the DM predicted 
by the thick disks in TC93 and M2. The difference between 
the model predictions at infinity increases at very small and 
very large Galactic longitudes, see Fig. [S] The reason for 
this difference is that while the exponential disk in model 
M2 has a uniform mid-plane electron density, the thick disk 
in TC93 shows a radial dependence of the electron density 
that is proportional to sech^ (radius) . 

4.5 Parameters for model M4 

In Fig.|6]we plot the DM at infinity, DMoo, against the EMoo 
that we calculated from observed Ha intensities fij |4.1|) . The 
DMoo were calculated as the sum of the prediction by model 
M2 at infinity plus the difference between the residual DMs 
and the model prediction by the thick disk in M2 for these 
residual DMs. If model M4 can accurately predict this dif- 
ference then it can outperform M2. (We used model M2 to 
predict DM at infinity since we will show in Sect. [5] that 
M2 predicts DM most accurately of all the models that we 
tested) Since we do not know a priori which ISM component 
is responsible for the the difference between the residual DM 



Figure 6. DMoo vs. EMqq for the lines of sight towards the se- 
lected pulsars from Table [T] (crosses) and globular clusters from 
Table [2] (circles) . The measured DMs were extrapolated to infin- 
ity using model M2 and include the remainder after fitting the 
thick disk in M2 to the residual DMs. The E Moo were calculated 
using the E(B-V) from ISchlegeT et al The best bisector 

fit DMoo = 22.67 EMj^"^* is indicated by a red line. 



and the DM predicted by the thick disk in M2, we cannot 
predict how this term changes when we extrapolate it to 
infinity. 

Interestingly, DMoo and EMoo show a correlation over 
2 orders of magnitude, for both isolated pulsars and pul- 
sars in globular clusters. The existence of this correlation 
is not a new result, since it was already described by 
iBerkhuiisen et a l. (2006). With an OLS bisector m ethod (see 
the instructive review in iFeigelsori fc Babul [l993 ) we found 
that the best fitting power law is described by 



DMoo = (22.67 ± 2.20) EM^*^"* =^ "'"^ (6) 

With our model M4 we investigate whether we can use this 
correlation to predict DM by using the observed EM as a 
proxy. 

We established the 1-cr errors of the coefficients of the 
DMoo-EMoo fit by running a Monte Carlo simulation of the 
distribution of the asymptotic values of DMoo.± • For each re- 
alisation we determined the scale height for the thick disk in 
M2 fij I4.4|l . which we then used to calculate DMoo from the 
measured DMs. By using a residual resampling bootstrap 
method we then determined the errors in the coefficients of 
the DMoo - EMoo fit, and we used the ensemble of all Monte 
Carlo realisations to establish the overall errors in these co- 
efficients (viz. the errors that include both the scatter in 
the asymptotic DMoo,±, and the scatter in the fitted DMoo 
- EMoo relation that is due to structure in the ISM). We 
checked that these errors did not change significantly when 
we changed the number of Monte Carlo realisations or the 
number of resamplings used in the bootstrap method. 



Modelling the Galactic distribution of free electrons 11 



5 COMPARISON BETWEEN MODEL 
PREDICTIONS AND THE DATA 

The predictions for DM made by the various models will not 
be perfect; instead they will scatter around the DMmodel = 
DMobserved line (Fig.[71 we will abbreviate 'DMobserved' to 
'DMobs' from here on). When analysing the scatter that the 
different models produce, one considers models with smaller 
fractional deviations to be better than models with small ab- 
solute deviations; we therefore plotted Fig. [7] with log axes. 
This figure also shows the ordinary least squares bisector 
fit of log(DMmodel) as a function of log(DMobs) as a red 
line. All models apart from TC93 clearly show a systematic 
deviation between the red line and the DMmodel — DMobs 
line (the dotted line). 

In this section we will investigate which model predicts 
the observed DMs most accurately over the entire range of 
observed DMs. Accuracy entails both the absence of system- 
atic deviations and a small scatter around the DMmodel = 
DMobs line. Although models may be very accurate over a 
small range in DM, where the red line intersects the dotted 
line, we will consider the predictive quality of the models 
over the full distance range, i.e., from nearby pulsars to pul- 
sars in distant globular clusters. 

One can define different ways for quantifying a (frac- 
tional) deviation between DMmodel and DMobs. An obvi- 
ous choice is DMmodel/DMobs, and other possibilities are 
(DMmodel-DMobs) /DMobs, and (DMmodel-DMobs)/(\/2 
DMobs), which respectively quantify the vertical distance 
and the perpendicular distance between a point and the 
line DMmodel=DMobs. As the latter two are both propor- 
tional to DMmodel/DMobs, the distribution that shows the 
smallest scatter in DMmodel/DMobs also shows the small- 
est scatter for the other 2 choices for quantifying deviations: 
the way one calculates the deviation between DMmodel and 
DMobs is not critical when defining what constitutes the 
best model. 

In Fig. [8]we show histograms of the distribution of DM- 
model/DMobs for the 45 lines of sight that were selected in 
§ 14.21 and we indicate data towards solitary pulsars with 
hatched bins, and data from pulsars in globular clusters 
with shaded bins. We used the subset of 45 lines of sight 
so that we can include model M4 in our comparison. Later 
in this section we will present results for all 68 available 
lines of sight. To emphasize the importance of fractional 
deviations between DMmodel and DMobs, we plotted the 
ratio DMmodel/DMobs using a logarithmic x-axis. System- 
atic deviations between the modelled and measured DMs 
show up in this figure as different distributions for the filled 
and hatched bins, since the pulsars in globular clusters have 
much larger DMobs than the solitary pulsars. The panels in 
this figure also show the normal distributions that were fit- 
ted to these data sets (red lines). The mean and spread of the 
fitted Gaussian curve were calculated from the distribution 
of DMmodel/DMobs of the 45 lines of sight; the amplitude 
of the fitted curve guarantees the correct normalization of 
the distribution. 

As is clear from comparing the histograms to the fitted 
normal distributions, the latter are not a very good descrip- 
tion of the data points. This means that the quality of the 
models can not be described in terms of a mean and spread 
of the fitted Gaussian. Also, the concept of statistical out- 



liers is ill-defined since sigma is not a good measure of the 
width of the distribution. Therefore, to compare the quality 
of the different models, we count for each model how many 
lines of sight have a DMmodel that deviates by less than a 
factor of N from DMobs (so N = DMmodel/DMobs when 
that ratio is larger than 1, otherwise N = DMobs/DMmodel) 
for different N - the result is the cumulative distributions of 
Figs.lQland lfOI where N is plotted along the x-axis. The best 
model should for each source have the smallest N compared 
to other models. This translates into a cumulative distribu- 
tion in Figs. [9]and[T0]that rises more steeply than the other 
distributions ('which lies above the other distributions' in 
these figures). This way of counting automatically incorpo- 
rates the systematic offsets between modelled and observed 
DMs that are seen in Fig. [7] Increasing the thick disk in 
TC93 has produced a model that predicts DM more accu- 
rately than the other models, both for the subset of 45 'good' 
sightlines, and for the full sample of 68 sightlines. 

In TableOwe list for each of the models the values of N 
within which 50%, 75% and 90% of the lines of sight lie. This 
table and Figs. [5] and [TO] show that all models can predict 
DM to within a factor of 1.5-2 for 3/4 lines of sight. The 
DMs of sightlines that lie further than 5° from the Galac- 
tic plane, and avoid local features in the Ha maps, can be 
predicted even more accurately, to within 40-70% for 3/4 
sightlines. M4 fails to be as accurate as the other models. 
The scatter around the DMoo(EMoo) relation probably pre- 
dicts a DMoo that, when corrected for the finite distance of 
the pulsar, produces a DM which deviates strongly from the 
observed DM. Model M4 might be improved upon, either by 
finding a second parameter that tightens the DMoo-EMoo 
fit, or by considering a different proxy. 

The difference in predictive quality between the differ- 
ent models becomes more apparent when considering the 
quality with which the DMs of 75% and 90% of the lines of 
sight are predicted. Excluding M4, the most accurate mod- 
els predict DM to within a factor of 2 of the observed DM, 
while the poorest performing models predict DM within a 
factor of 3 for 9/10 of sightlines. The DMs of the 23 sightlines 
that lie within 5° of the Galactic plane or that lie towards 
structures in Ha intensity are not by definition difficult to 
predict. It is somewhat surprising to find that the simpler 
models perform almost as well as the more complex models 
that are based on TC93 or NE200I. This is probably related 
to how the pulsars that we use to test the different models 
are distributed through the Milky Way. 



6 SIGHTLINES WITH ANOMALOUS N 

The cumulative distributions of N also lend themselves for 
defining the number of 'outliers' for which the modelled DMs 
deviate strongly from the measured DMs, for example by 
looking up the fraction of lines of sight with N > 2. In Table 
Owe list for each of the models which lines of sight have N 
that lie above the 90th percentile, supplemented with the 
lines of sight that have N > 2. Although in general differ- 
ent models have difficulty in predicting the DMs of differ- 
ent sightlines, most models cannot accurately predict the 
DM of J0358-h5413, J0630-2834, J0659-hl414, J0826-h2637, 
J0835-4510, JI643-1244, and J1909-3744. Only in the case 
of J 1909-3744 is the modelled DM larger than the observed 
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Figure 7. Comparison between the modelled and measured DMs (vertical and horizontal axes) for the 12 models from Table[3] Isolated 
pulsars are indicated by crosses, pulsars in globular clusters as open circles. The red line shows the OLS bisector fit to the data points. 
For model M4 the least-squares fitting algorithm did not converge (which is probably related to the very large scatter of the data points) 
so that no OLS bisector line was drawn in that case. 
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Figure 8. Histograms of the ratios between tlie modelled and observed DMs for each of the 12 models from Table [3] The hatched 
and shaded bins indicate data from isolated pulsars and from pulsars in globular clusters respectively. The red curve shows the normal 
distribution whose mean and width were determined using robust statistics, the amplitude of the curve then follows from a least-squares 
minimisation. 
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Figure 9. Cumulative distribution of the number of linos of sight 
where the modelled DM lies within a factor of N of the observed 
DM. We only plot N for the 45 lines of sight selected in ^ l4.2l so 
that model M4 can be included, and we used 3 panels to avoid 
cluttering the plot. 

DM, for all models. The DM that were calculated for the 
other sightlines are smaller than the observed DM. 

When the high-resolution Ho? maps from SHASSA or 
VTSS indicate an Hll region in the direction of one of our 
lines of sight, we can use that line of sight to put a upper 
or lower limit on the distance to the Hll region. If the ob- 
served DM is larger than the predicted DM, the real ISM 
contains more structure than is included in the model, which 
puts the pulsar or globular cluster beyond the Hll region. 
Similarly, when the sightline does not show an anomalous 
N, the pulsar or globular cluster must lie in front of the 
H 11 region. This way the distances that were determined for 
J0358-H5413, J0659-H414 and J1643-1244 give upper hmits 
to the H 11 regions that must lie in front of them. The prox- 
imity of J0358-f 5413 to the Galactic plane, combined with a 
sightline of over 1 kpc, could mean that more than one H ll 
region contributes to its DM. Although J0835-4510 is seen 
towards the Gum and Vela nebul ae, its distance places it 
in fro nt of these nebulae (see also iMitra fc RamachandranI 
I2OOII ). It turns out that models which are based on NE2001 
unrferpredict the DM of this pulsar, while the other models 
overpredict the DM. 
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Figure 10. The equivalent of Fig. [9] for all 68 available lines of 
sight from Tables [T] and |2] 

We analysed the high-resolution SHASSA and VTSS 
maps of Hq intensity towards J0630-2834, but we did not 
find strong Ha emission in the vicinity of these lines of sight. 
J0826-I-2637 is only covered by the low-resolution WHAM 
data. There is however diffuse Hq emission over a wider 
region surrounding the lines of sight towards J0630-2834 and 
J0826-I-2637, but to identify such regions would require more 
advanced selection tools. 

We did not exclude these lines of sight from our analysis 
since we want to keep our selection criteria simple. In prac- 
tice, since these lines of sight have the highest N in most 
models, they affect the cumulative distributions of these 
models in the same way, and only above the 90th percentile. 



7 CONCLUSIONS 

In this paper we have analysed how well twelve different 
models of the Galactic free electron density rie predict the 
measured DMs of pulsars at well-known distances. These 
distances were derived from either parallax measurements, 
or from association of the pulsar with a globular cluster. 
Eight of these models have already been published in the 
literature. We determined a scale height of 1.6 kpc for our 
exponential thick disk model Ml, where we used only sight- 
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Table 5. Overview of sightlines that lie above the 90th percentile of the cumulative distribution of N (values of N indicated in italics), 
supplemented with sightlines for which N > 2 (i.e. DMjnodel > 2xDM_obs, or DMjnodel < l/2xDM_obs). Name : line of sight was not 
part of the 45 lines of sight that were selected in § 14.21 Since we are investigating the N of all 68 lines of sight we left out the results 
for model M4, which can only predict DMs well for the subset of 45 lines that were selected in the aforementioned section. In the final 
column we list the Galactic latitudes of the lines of sight, rounded off to the nearest integer. 
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lines with |6| > 32.5° to avoid the more complex structure 
at smaller Galactic latitudes. In models M2 and M3 we re- 
placed the t hick disk from t h e mod els bv lTavlor fc Corded 
l|l993l ') and ICordes fc Lazid l|2002l ) by exponential thick 
disks. First we had to correctly account for the DM contribu- 
tions by the other components in these models, which we did 
by subtracting these DMs from the observed DMs before fit- 
ting a new scale height for the thick disk. The thick disks in 
M2 and MS have scale heights of 1.6 and 1.3 kpc, resp. When 
we limit the radial extent of the thick disk in model M2 to 
27.7 kpc, it produces for !&!=0° the same DM at infinity as 
the or iginal thick disk from the model bv lTavlor fc Corded 
Although this radius on average reproduces the DM 
of the original thick disk, the DM difference between the old 
and new thick disk varies somewhat as a function of Galac- 
tic longitude. With model M4 we explore using the emission 
measure as a proxy to predict the observed scatter in DM 
that the other models cannot explain. 

A few of the components in the model by 



ICordes fc Lariol (|2002l ') contribute also to the DMs of 
sightli nes that lie furt h er tha n 40° from the Galactic plane. 
Since ICaensler et al.l (|2008l ) only fitted an exponential 
thick disk, incor porating their model for the thick disk 
in the model by ICordes fc Lazid l|2002l ) leads to internal 
inconsistencies. 

We developed a novel approach for comparing the pre- 
dictive quality of the different models, by first calculating for 
each sightline the factor N by which the modelled DM devi- 
ates from the observed DM (so N = modelled DM/observed 
DM when that ratio is larger than 1, otherwise N = observed 
DM/modelled DM), and then comparing the cumulative dis- 
tributions of N between the models. Most models can predict 
DM to within a factor of 1.5-2 of the measured value fo r 3/4 
lines of sig ht, and the model bv lTavlor fc Corded (|l993l ) that 
we updated with a more extended thick disk consistently 
performs better than the other models that we tested. Our 
fourth model predicts DM least accurately. Surprisingly, the 
simpler models that we tested perform almost as well as the 
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more complex models. This is probably related to the sparse 
spatial distribution of the sightlines that we used. Finally, 
by comparing which pulsar DMs are difficult to predict by 
most models we identified pulsars that lie beyond Hii re- 
gions. The accurate pulsar distances then place upper limits 
on the distances to these H ii regions. 

The analysis we present here will benefit from having 
a larger sample of sightlines with known distances and dis- 
persion measures, which can come from either parallax mea- 
surements of additional pulsars, or from discoveries of more 
globular clusters with pulsars in them. A different approach 
in the modelling philosophy, for example by including Ha 
data or a different proxy, warrants further investigation, as 
it can help explain the detailed structure in DM. 
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